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ABSTRACT 

We present Hi observations of the Sculptor Group starburst spiral galaxy NGC 253, ob¬ 
tained with the Karoo Array Telescope (KAT-7). KAT-7 is a pathfinder for the SKA pre¬ 
cursor MeerKAT, under construction. The short baselines and low system temperature of the 
telescope make it very sensitive to large scale, low surface brightness emission. The KAT-7 
observations detected 33% more flux than previous VLA observations, mainly in the outer 
parts and in the halo for a total Hi mass of 2.1 ± 0.1 x 10® M©. Hi can be found at large 
distances perpendicular to the plane out to projected distances of ^9-10 kpc away from the 
nucleus and ~13-14 kpc at the edge of the disk. A novel technique, based on interactive pro¬ 
file fitting, was used to separate the main disk gas from the anomalous (halo) gas. The rotation 
curve (RC) derived for the HI disk confirms that it is declining in the outer parts, as seen in 
previous optical Fabry-Perot measurements. As for the anomalous component, its RC has a 
very shallow gradient in the inner parts and turns over at the same radius as the disk, kinemat¬ 
ically lagging by ~100 km/sec. The kinematics of the observed extra planar gas is compatible 
with an outflow due to the central starburst and galactic fountains in the outer parts. However, 
the gas kinematics shows no evidence for inflow. Analysis of the near-IR WISE data, shows 
clearly that the star formation rate (SFR) is compatible with the starburst nature of NGC 253. 

Key words: techniques: interferometric - galaxies: individual: NGC 253 - galaxies: starburst 
- galaxies: kinematics and dynamics - galaxies: haloes 


1 INTRODUCTION 

The seven-dish KAT-7 array (Carignan et al. 2013) was built as an 
engineering testbed for the 64-dish Karoo Array Telescope, known 
as MeerKAT, which is the South African precursor of the Square 
Kilometre Array (SKA). KAT-7 and MeerKAT are located close 
to the South African SKA core site in the Northern Cape’s Karoo 
desert region. Construction of the array was completed in Decem¬ 
ber 2010. The array is extremely compact, with baselines ranging 
from 26 m to 185 m and the receivers have a very low Tsys ~ 30K. 
While its main purpose is to test technical solutions for MeerKAT 
and the SKA, scientific targets, such as NGC 253, were also ob¬ 
served during commissioning to test the Hi line mode. In this pa¬ 
per, we present over 150 hours of observations taken with KAT-7 
in order to study its large scale extra-planar Hi gas component. 

NGC 253 is the nearest example of a galaxy undergoing en¬ 
hanced star formation in the South, M82 being the nearest one in 
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the North. Its nucleus contains a starburst with a star-formation rate 
(SFR) of ~5 Mq yr“^, roughly 70% of the rate of the entire galaxy 
(Wik et al. 2014). It is considered a prototype starburst galaxy (see 
e.g. Sakamoto et al. 2011; Strickland et al. 2002) and Rieke, Lebof- 
sky & Walker (1988) suggest that it may be in an earlier starburst 
phase than M82. Of particular interest is that, despite the high SFR, 
NGC 253 is a non-interacting (non-merger) system, which does not 
exclude that it occurred in the past (see Sec. 7.4). 

One very interesting feature is the diffuse extra-planar Hi sur¬ 
rounding the X-ray emission, first seen in the ATCA observations 
of Boomsma et al. (2005). Not only is the nuclear starburst form¬ 
ing stars at a fairly high rate but it is also thought to produce a 
super-wind (Heckman, Armus Sc Miley 1990), which expels ma¬ 
terial into the halo. The high sensitivity of KAT-7 to large scale 
diffuse emission (Carignan et al. 2013) allows us to better study 
the known extra-planar gas component. Such extra-planar gas is 
expected when looking at a deep UV image showing emission far 
from the plane, compared to a 2MASS image (Jarrett et al. 2003), 
which shows the main stellar disk component (Figure 1). 
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Figure 1. 2MASS K-band (top) and GALEX FUV (bottom) images of NGC 
253. 


In NGC 253, most of the activity is confined within ~1 kpc of 
the nucleus, fuelled by bar-driven dense cold molecular gas that ob¬ 
scures the starburst (Rieke et al. 1980; Hougton et al. 1997; Engel- 
bracht et al. 1998; Sakamoto et al. 2011). Recent X-ray, optical and 
ALMA mm-wave observations trace outflowing gas that is powered 
by the starburst (Mitsuishi, Yamasaki & Takei 2013; Bolatto et al. 
2013). Owing to proximity, the two nearest starbursts, NGC 253 
and M82, are ideal laboratories for detailed (pc-scale) study of nu¬ 
clear outflows and galactic-disk fountains. Conversely, their global 
properties - which place the starbursts in context with the total gas 
and stellar content - are a challenge to measure due to their large 
angular size and brightness (frequently saturating detectors). While 
large angular scales are not a problem at optical wavelengths, for 
radio interferometry imaging studies, one needs the proper short 
baselines to be able to detect large structures and the compactness 
of KAT-7 makes it better suited than e.g. the VLA for the observa¬ 
tion of those large scales in nearby objects (see Table 2). 

Recent deep Hi observations of several nearby galaxies indi¬ 
cate that up to 15% of the neutral hydrogen of a spiral galaxy is 
located in the halo (Heald et al. 2014). The gas outside the main Hi 
disk of spirals has been studied in many galaxies: NGC 891 (Swa- 
ters, Sancisi & van der Hulst 1997; Oosterloo, Fratemali & Sancisi 
2007), NGC 2403 (Schaap, Sancisi & Swaters 2000; Fratemali et 
al. 2001), NGC 2997 (Hess et al. 2009), NGC 3198 (Gentile et al. 
2013), NGC 4244 (Zschaechner et al. 2011), NGC 4559 (Barbi- 


eri et al. 2005), NGC 4565 (Zschaechner et al. 2012), NGC 5775 
(Lee et al. 2001), NGC 6946 (Boomsma et al. 2008), M31 & M33 
(Westmeier, Braun & Thilker 2005; Thilker et al. 2004) and the 
Milky Way (MW) through High Velocity Clouds (HVCs) (Wakker 
& van Woerden 1997), the Magellanic Stream (Mathewson, Cleary 
& Murray 1974) and the Leading Arm (Putman et al. 1998). See 
also Heald et al. (2011, 2014) for a description of the HALOGAS 
survey and Sancisi et al. (2008); Putman, Peek & Joung (2012) for 
reviews. 

The interaction between the gas in the halo and the disk of spi¬ 
rals is believed to play an important role in their evolution. Halo gas 
connects the baryon-rich intergalactic medium (IGM) to the star¬ 
forming disks of galaxies (Putman, Peek & Joung 2012). The clas¬ 
sical "galactic fountain" scenario (Schapiro & Field 1976), where 
hot gas is expelled from the disk through "galactic chimneys" cre¬ 
ated by multiple supernova explosions from clusters of young mas¬ 
sive stars (Normandeau, Taylor & Dewdney 1996) can explain the 
presence of hot gas in the halo. It is believed that, as the gas ex¬ 
pands, it eventually converts back to Hi as it cools through radiative 
losses; raining back down on the disk (Bregman 1980) and feed¬ 
ing subsequent star formation (see also the models for the MW: 
Marasco, Fratemali & Binney 2012). A good example of such ex¬ 
changes between the disk and the halo is the Ha kinematical study 
by Cecil et al. (2001) of one of the filament in the super wind out¬ 
flow of NGC3079. Long ago, Oort (1966) also suggested an exter¬ 
nal origin for some of that gas in the form of primordial gas clouds 
left over from the formation of the galaxies. Finally, as suggested 
by the Magellanic Stream and the Leading Arm in the MW, some 
of that halo gas could come from ISM torn out of dwarf galaxies 
during close encounters with a large spiral. 

The optical and IR parameters of NGC 253 are summarized 
in Table 1. It is the brightest and earliest type among the five 
late-type spirals composing the nearest group to us, the Sculptor 
group (Puche & Carignan 1988). It is classified SAB(s)c by de Vau- 
couleurs et al. (1991) but, as can be seen in Figure 2, it is clearly 
a barred SB(r)c galaxy on IR images. The different recent distance 
estimates from the Planetary Nebulae Luminosity Function (PNLF) 
and the Tip of the Red Giant Branch (TRGB) distance indicators 
(Rekola et al. 2005; Mouchine et al. 2005; Dalcanton et al 2009; 
Jacobs et al. 2009; Hislop et al. 2011) give a distance of 3.5 Mpc, 
for a scale of ~1 kpc per arcmin. While NGC 55 and NGC 300 
are on the near side and NGC 7793 on the far side of the Sculp¬ 
tor group, NGC 253 is at the mean distance of the group with NGC 
247, which is only at a projected distance of ~350 kpc to the North. 

The first detailed Hi observations of this galaxy (Huchtmeier 
1972; Combes, Gottesman & Weliachew 1977) show clearly an 
asymmetric Hi distribution. The presence of a strong continuum 
source in the central regions (Hummel, Smith & van der Hulst 
1984) suggests that any kinematical analysis in the very inner parts 
is to be undertaken with great care. From the low spatial resolution 
single dish early Hi observations, the rotation curve was believed to 
be declining at large radii (Huchtmeier 1975) but higher resolution 
VLA observations (Puche, Carignan & van Gorkom 1991) show a 
rotation curve still rising at the last observed velocity point (~12'). 
However, recent Fabry-Perot observations (Bland-Hawthom, Free¬ 
man & Quinn 1997; Hlavacek-Larrondo et al. 2011) detected Ha 
emission further out than previous Hi observations. While the Ha 
kinematics agree with the Hi for R < 12', it suggests a declining ro¬ 
tation curve between 15'and 19'. It will be interesting to see if our 
KAT-7 observations extend to those radii and confirm this decline. 

Figure 3 shows previous higher spatial resolution (68" x 68" 
vs 213" X 188") VLA observations (Puche, Carignan & van 
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5 arcmin 


Figure 2. Wise W1 (blue), W2 (green), W3 (orange) and W4 (red) color- 
combination image ofNGC 253. 


Parameter 

Value 

Morphological type^ 

SB(r)c 

Right Ascension^ (J2000) 

0^* 47“ 33.L 

Declination^ (J2000) 

00 

0 

IT) 

Distance (Mpc) 

3.47 ± 0.14 

Scale (kpc arcmin~^) 

~ 1 

Isophotal major diameter^, D 25 

27.5' 

Holmberg radius^, Rho 

17.5' 

Axis ratio^ (q=b/a), R 25 

0.25 

Inclination (qo=0.15) 

78.3° 

Position Angle 

52° 

Total apparent B magnitude^ 

8.04 

Corrected apparent B magnitude^ 

7.02 

Absolute B magnitude 

-20.68 

Total B luminosity, Lq 

2.78 X 10°° 

Absolute W 1 magnitude, WISE 

-24.24 

Color (W1 - W2), WISE 

0.21 

Log(M.), WISE, Mq 

10.33 


(1) See Figure 2 

(2) de Vaucouleurs et al. (1991) 

(3) Radius at fj,B = 26.6 mag arcsec 


Table 1: Optical & IR (WISE) parameters ofNGC 253. 

Gorkom 1991). Note however that the final KAT-7 data have a 
factor of ~4 better velocity resolution. Regarding the Hi distribu¬ 
tion, the first feature to notice is the strong absorption in the centre, 
which prevents any reliable kinematics to be extracted close to the 
starburst region. It can be seen that besides tracing the spiral arms, 
the Hi is really concentrated around the nucleus. As for the veloc¬ 
ity field, it is a textbook case in terms of regularity. Only a slight 
increase in position angle is suggested in the very outer parts. The 
most interesting map is the velocity dispersion map, which shows a 
clear increase ~45-50 km s“^ around the starburst region while it 
is fairly constant ~10-15 km s“^ in the rest of the disk. Because of 
the short exposure time (~30 minutes in D configuration) and the 
lack of short spacings, those observations only show the disk emis¬ 


(1) 

Telescope 

(2) 

Min. bl 
(m) 

(3) 

Ant. diam. 
(m) 

(4) 

LAS 
(min. bl) 

(5) 

LAS 

(sha. lim.) 

KAT-7 

26 

12 

16.7' 

36.1' 

VLA-D 

35 

25 

12.4' 

17.3' 

WSRT-MS 

36 

25 

12.0' 

17.3' 

GMRT 

100 

45 

4.3' 

9.6' 

ATCA 

31 

22 

14.0' 

19.7' 


Table 2: Largest Angular Sizes (LAS) of different synthesis tele¬ 
scopes. Given for the different telescopes (1) are: minimum base¬ 
line (2), antenna diameter (3), LAS from the minimum physical 
baseline (4), and LAS down to the shadowing limit (5). 


sion and were insensitive to the large scale (> 15') low column 
density extra planar Hi. 

In theory, the Largest Angular Sizes (LAS) accessible to an 
aperture synthesis telescope is given by: 

LAS = 0.6 X A/Binin (1) 

When trying to observe LAS, Table 2 illustrates why KAT-7 
should be a good instrument for this purpose. If we consider the 
shortest physical baseline, it can be seen that the gain of KAT-7 is 
~33% over the VLA, WSRT (in its Maxishort configuration) and 
ATCA. In fact, the minimum projected baseline is the antenna di¬ 
ameter. This means that if we instead consider the shadowing limit, 
the gain over the VLA, WSRT and ATCA is a factor of 2 (~ 36' vs 
18'). 

In practice, the largest scale that can be reliably imaged cannot 
be assessed as simply as Table 2 suggests. When less flux is recov¬ 
ered with an interferometer than with a single dish, the shortcoming 
is often explained by a lack of short spacings. Yet it may often be 
the case that additional care in calibration and imaging can resolve 
the discrepancy. The resulting quality when imaging large-scale 
emission is critically dependent on several factors including the im¬ 
age reconstruction techniques that are employed (classic CLEAN, 
MEM, multi-resolution CLEAN, etc). It is important to note that 
the spatial scales of interest should be assessed per spectral chan¬ 
nel, not for the entire galaxy. In single channels, the angular scales 
are always more restricted. A good illustration of imaging struc¬ 
tures larger than the nominal LAS is the study of NGC 628 by 
Kamphuis & Briggs (1992), in which features as large as ~25'are 
seen in individual channels (see their Lig. 1 showing the naturally 
weighted channel maps), although emission on scales larger than 
~17.5'should not be seen by the VLA in D configuration. Nev¬ 
ertheless, it remains true that KAT-7’s smaller dish diameter and 
shorter baselines enable it to image larger structures than with the 
VLA, WSRT and ATCA. 

The remainder of this paper is as follow. In Sec. 2, a descrip¬ 
tion of the KAT-7 observations and data reduction will be given. 
Sec. 3 will discuss the continuum emission and Sec. 4 will describe 
the Hi content and distribution. Sec. 5 will explain how to distin¬ 
guish between the Hi in the disk and the Hi in the halo while Sec. 6 
will study the Hi kinematics of both components. The main results 
will be discussed in Sec. 7 and, finally, a summary and final con¬ 
clusions will be given in Sec. 8. 
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Figure 3. VLA moment maps (top: Hi distribution; middle: velocity field; 
bottom: velocity dispersion) of NGC 253 derived from the data set pre¬ 
sented by (Puche, Carignan & van Gorkom 1991). The velocity resolution 
is 20.6 km s“^ and the spatial resolution is 68” X 68”. 


Parameter 


Value 

Number of antennas 


7 

Dish diameter 


12 m 

Min baseline 


26 m 

Max baseline 


185 m 

Frequency range 

1200 - 1950 MHz 

Max instantaneous bandwidth 

256 MHz 

Polarisation 


Linear H & V 

T 

-isys 


26 K 

Aperture efficiency 


0.65 

System Equivalent Flux Density 

1000 Jy 

Latitude 


-30:43:17.34 

Longitude 


21:24:38.46 

Elevation 


1038 m 

Digital back-end 

ROACH boards 

Table 3: KAT-7 specifications. 

mode total BW 

number of 

channel BW 

(MHz) 

channels 

(kHz) 

cl6n2M4k 1.5625 

4096 

0.381 

cl6n7M4k 6.25 

4096 

1.526 

cl6n25M4k 25 

4096 

6.104 

cl6n400M4k 256 

1024 

390.625 


Table 4: KAT-7 correlator modes. 


2 KAT-7 OBSERVATIONS AND DATA REDUCTION 

The KAT-7 dishes have a prime-focus alt-az design with a F/D of 
0.38, optimized for single-pixel L-band feeds. The low noise am¬ 
plifiers (LNAs) for the feeds are cryogenically cooled to 80 K using 
Stirling coolers. The key system specifications for KAT-7 are sum¬ 
marized in Table 3. The digital backend of the system uses the Re- 
configurable Open Architecture Computing Hardware (ROACH- 
https://casper.berkeley.edu/wiki/ROACH), which is a flexible and 
scalable system enabling spectral line modes covering a wide range 
of resolutions. Table 4 gives the details of the recently commis¬ 
sioned correlator modes. Digital filters give a flat bandpass over the 
inner 75% of the band with a rapid roll-off at the edges of the band. 

The parameters of the KAT-7 observations of NGC 253 are 
given in Table 5. The data were collected over 17 observing ses¬ 
sions between 2013 March and 2013 July using the cl6n25M4K 
spectral line mode. This correlator mode gives velocity channels 
of 1.3 km s“^ over a flat bandpass of 4000 km s“^. The larger 
bandwidth was used to look for background galaxies in the field 
and to ensure a good estimation of the continuum for subtraction. 
All antennas were in the array during all the observing sessions. 
The median time of each session on target was 9 hours for a total 
project time of 153.5 hours, including calibration and slew time. 
Total time on NGC 253 was 115.6 hours. 

Each of the seventeen observing sessions were reduced sep¬ 
arately. All data calibration was done using standard calibration 
tasks in the Common Astronomy Software Applications (CASA 
3.4.0) package (McMullin et al. 2007). The KAT-7 primary beam 
is large enough to observe NGC 253 in a single pointing. Phase 
drifts as a function of time were corrected by means of a nearby 
point source (0023-263) observed every ten minutes. This source 
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Parameter 

Value 

Start of observations 

5 March 2013 

End of observations 

19 July 2013 

Total integration (on source) 

115.6 hours 

FWHM of primary beam 

1.08° 

Total Bandwidth 

25 MHz 

Central frequency 

1419 MHz 

Channel Bandwidth 

6.1 kHz 

Number of channels 

4096 

Channel width 

1.28 km s“^ 

Maps gridding 

20" X 20" 

Maps size 

257 X 257 

Flux calibrator 

3C138 

Phase/bandpass calibrator 

0023-263 

Robust = 0 weighting function 
FWHM of synthesized beam 
RMS noise 

Column density limit (3(7) 

213"x 188" 
1.0 mJy/beam 
1.3 X 10^® cm"2 


Table 5: Parameters of the KAT-7 observations. 


was also used to correct for variations in the gain as a function of 
frequency (bandpass calibration). The absolute flux scale was set 
by observations of 3C138. Comparisons of the flux measurements 
on the observed calibrators suggests that the absolute flux uncer¬ 
tainties are on the order of 5%. Variations in the bandpass are on 
the order of 1%. 

Initial imaging revealed which channels were free of HI emis¬ 
sion. Continuum emission was subtracted from the raw UV data 
by making first order fits to the line free channels using the CASA 
task UVCONTSUB. The calibration was then applied and NGC 
253 was then SPLIT from the calibration sources. In order to make 
sure no residual calibration errors remained in the data, a test cube 
was made using the CASA task CLEAN. KAT-7 does not employ 
doppler tracking and CASA does not fully recognize frequency 
keywords, so special care was taken to produce uv data sets and 
test cubes with the proper velocity coordinates (see Carignan et al. 
2013). The individual calibrated continuum subtracted uv data sets 
were then combined together using the CASA task CONCAT. 

Preliminary imaging of the combined data in CASA revealed 
the presence of artifacts in the form of horizontal lines which we 
were able to identify as the accumulation of RFI at u=0. This is 
most likely due to the fact that for those visibilities the fringe rota¬ 
tion is zero. These were removed by flagging all visibilities near u = 
0. This uncovered further artifacts in the form of diagonal lines par¬ 
allel to the major axis of NGC 253. To remove this problem, it was 
found necessary to use the task SELFCAL in Miriad (Sault, Teuben 
& Wright 1995), on the strong continuum of NGC 253. This re¬ 
moved most of these artifacts giving a final rms of 1 mjy/beam 
with channels smoothed to 5 km s“^ in the final cube. 


3 CONTINUUM EMISSION 

Figure 4 shows the KAT-7 continuum map of NGC 253. The to¬ 
tal flux at 1.419 GHz is 6.8 ± 0.3 Jy. This is comparable with the 
previous measurements of 6.5 ± 0.6 Jy at 1.415 GHz (de la Beau- 
jardiere et al. 1968) and 6.7 ± 0.7 Jy at 1.430 GHz (Heidmann 
1963). When comparing Fig. 4 to the observations at 1.46 GHz of 


Hummel, Smith & van der Hulst (1984) or at 1.49 GHz of Con¬ 
don (1987), the emission is slightly more extended (~20%) along 
the major axis but 60% more along the minor axis (see also Car- 
illi et al. 1992). In fact, the emission reaches the very edge of the 
optical disc on this very deep image from Malin & Hadley (1997). 
However, the difference between our continuum map and the one of 
Hummel, Smith & van der Hulst (1984) is not surprising since they 
themselves mentioned that their emission is sitting in the middle of 
a negative bowl, a signature of the lack of short enough baselines 
at the VLA. 

Looking at the very high flux in the centre, it is clear that the 
starburst is responsible for a large part of the emission, while the 
spiral arms are more responsible for the lower flux in the outer disk. 
Along the major axis, it can be seen that the emission stops at the 
end of the spiral arms. Contrary to what is seen in the observations 
of some edge-on galaxies, where the axis ratio of the continuum 
emission is similar to the axis ratio of the optical disk (Hummel, 
Sancisi & Ekers 1984), the axis ratio is twice that of the optical 
disk (0.50 vs 0.25) in the case of NGC 253. 

However, NGC 253 is not the only exception. NGC 4631 (Ek¬ 
ers &. Sancisi 1977), which has strong X-ray emission originating 
from the nucleus and is interacting with a dwarf companion along 
the minor axis, also has a much thicker radio disk. This is only 
beaten by NGC 253’s twin M82 (Seaquist & Odegard 1991), which 
has an axis ratio ~L0 (vs optical ~0.4). Higher resolution (VLA, 
A configuration) 20cm observations of the centre of NGC 253 (Ul- 
vestad & Antonucci 1997) show an even thicker disk with an axis 
ratio ~0.65. 

Allen, Baldwin & Sancisi (1978) explain that the sources of 
relativistic electrons are distributed like the light in the optical disk 
and that the electrons diffuse outwards. However, the larger thick¬ 
ness of the radio disk in this case suggests that the starburst may be 
responsible for an extra z-component of the magnetic field, which 
enhances the propagation of the relativistic electrons in the direc¬ 
tion perpendicular to the plane. As discussed by Heald (2012), 
a substantial vertical displacement of an initially plane-parallel- 
ordered magnetic held may be driven by a star formation event. A 
detailed study of the magnetic field in the nuclear outflow of NGC 
253 can be found in Heesen et al. (2011). 

Since the KAT-7 observations were carried out over several 
months, we have a unique opportunity to investigate whether or not 
the continuum in the central regions exhibits time variability. To 
this end we calculate the continuum flux inside a 2' radius cover¬ 
ing the absorption feature for all seventeen separate data sets. No 
significant variation in continuum flux is observed. All measured 
fluxes are within the measurement uncertainties (~ 5 %). It should 
be noted here again that the bulk of the radio continuum comes 
from the central starburst, reducing the likelihood of detecting sig¬ 
nificant variations in continuum flux. 


4 Hi CONTENT AND DISTRIBUTION 

The total Hi distribution map, shown in Figure 5, was derived using 
the task MOMNT (moment 0) in AIPS (Greisen 2003). It is super¬ 
posed on a composite image of FUV from Galex, IR from WISE 
and soft X-rays from ROSAT (0.1-0.4 keV). The faintest level goes 
down to ~ 1.0 X 10^® cm~^. At that level, the galaxy has a diame¬ 
ter of 34 ±2 kpc (the error being defined by the beam size), which 
is comparable to the optical diameter (cf. Tab. 1). It can be seen that 
a significant fraction of the Hi can be found away from the plane 
out to projected distances of ~9-10 kpc in the centre and 13-14 kpc 
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Figure 4. Continuum map of NGC 253, overlaid on a deep optical image 
(Malin & Hadley 1997). The contours are 7.5 (3 (t), 15, 30, 60, 120, 240, 
480, 960 mJy/beam. The beam is in the bottom left corner. 



Figure 6. Global profile (black) for the KAT-7 Hi observations with the 
kinematical Vsys indicated (see Sec. 6). Also plotted are the Hi profiles for 
the receding (red dotted) and approaching (blue dashed) components of the 
Halo (see Sec. 5). The channel width is 5 km s“^ . 


at the edge of the disk, which is comparable to the extent of the 
extra planar emission seen in the edge-on spiral NGC 891 (Ooster- 
loo, Fratemali & Sancisi 2007). Looking at Fig. 5, the FlI appears 
to surround the hot gas component in the outer parts. 

The global Fll profile of NGC 253 is given in Figure 6 . This 
was obtained using the task BLSUM in AIPS after primary beam 
correction. The asymmetry is clear with more Hi on the receding 
(SW) than on the approaching side. From it, a mid-point velocity 
of 242±4 km s“^ (50% level) is found. Because of the asymmetry, 
this value has to be preferred to the intensity weighted mean ve¬ 
locity as indicative of the systemic velocity. This can be compared 
to 243±2 km s“^ found by the HIPASS survey (Koribalski et al 
2004). However, for the kinematical analysis of Sec. 6 , we will use 
the kinematically determined value. The profile widths at the 20% 
and 50% levels are AV20 = 439± 4 km s“^ and AV50 = 412± 4 
km s“^, slightly larger than the Koribalski et al (2004) values. A 
total Hi flux of 728 ± 36 Jy.km s“^ is found (after primary beam 
correction), which corresponds to an Hi mass of 2.1 ± 0.1 x 10® 
Mq for a Mhi/Lb — 0.1 at our adopted distance of 3.5 Mpc. This 
is similar to the mass of 2.0 x 10® Mq found by Boomsma et al. 
(2005). Note that in the next section, we will separate the disk and 
the kinematically anomalous (halo) components. The mass of the 
anomalous gas is given in Table 6 . 

The radial Hi profile of NGC 253 is given in Figure 7 and 
compared to the profile derived by the previous VLA observations. 
It can be seen that, while the VLA was detecting the bright Hi disk 
(down to 2.4 x 10®® cm“®), KAT-7 detects ~33% more flux (down 
to 1.3 X 10^® cm“®), mainly in the outer disk and in the halo. This 
is not only coming from the longer integration but also because 
KAT-7 can detect scales larger than 15^ scales invisible to the VLA 
because of the lack of short spacings. 

The derived Hi parameters are summarized in Table 6 . While 
the Hi disk of a late-type spiral is often much more extended than 
the stellar disk (D 25 ), this is not the case for NGC 253. At a level 
of 1 Mq pc“® (1.25 X 10®° cm"®), the Hi disk of NGC 253 is 
equivalent to the stellar disk (isophotal major diameter, D25) and 
even at a level of 10^® cm"®, it is only equivalent to the Holmberg 
diameter (Dho ’■ diameter at a surface brightness level of /rs = 26.6 
mag arcsec"®). For example, NGC 300, another late-type spiral 
in the Sculptor group, has an Hi disk 50% larger than the optical 
disk (Puche, Carignan & Bosma 1990). What is exceptional with 



Figure 7. Comparison of the VLA (Puche, Carignan & van Gorkom 1991) 
and the KAT-7 Hi radial profiles. 

NGC 253 is not its radial extent but more its extent perpendicular 
to the disk. Despite the apparent truncation of the Hi disk, NGC 
253 is not Hi deficient with logfMni) = 9.32, slightly larger than 
<log(MHi)> for an isolated field Sc of 9.24 (Solanes, Giovanelli 
& Haynes 1996). 


5 SEPARATING THE DISK AND HALO Hi GAS 

The Hi total intensity map of NGC 253 clearly shows some of the 
Hi to be spatially separated from the main disk. Furthermore, the 
kinematics of this spatially extended Hi component differ from that 
of the main disk. Position-Velocity (PV) slices, extracted parallel to 
the major axis of NGC 253, provide clear evidence of kinematically 
anomalous Hi. Such a PV slice, passing through the centre of the 
galaxy, is shown in Figure 8 . For the approaching (NE) half of the 
galaxy, it clearly highlights the presence of an anomalous Hi com¬ 
ponent spanning a range of radial velocities that is offset from the 
high-intensity emission of the main disk; toward the systemic ve¬ 
locity of the system. Such anomalous Hi features are commonly 
seen in deep observations of nearby galaxies: e.g. NGC 891 (Oost- 
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Figure 5. Combined WISE (RGB) + FUV (magenta) + X-rays (yellow: Rosat) -H HI contours of the moment 0 map (scale: l'~ 1 kpc). The contours correspond 
to 1.3, 2.6, 5.2, 10.4, 20.8, 41.6 ,83.2, 166.4 & 253.3 xlO^^cm-^. The ROSAT data are from Pietsch et al. (2000). 


Parameter 

Value 

Mhi (total) 

2.1 ±0.1 xlO® M 0 

Mjqj (anom.)^ 

7.8 X 10'^ M 0 

Mpfj (anom.)/Mpff (total)^ 

3.5% 

AVao, ( 20 % level) 

439 ±4kms"^ 

AVso, (50% level) 

412 ±4 km s"^ 

Vsya (mid-point 50%) 

242 ±4 km s"^ 

Vsys (kinematical, TR) 

238 ±4 km s"^ 

Dao, (diam. at 10 ^° cm“^ level) 

29 ±2 kpc 

Dig, (diam. at 10^® cm“^ level) 

34 ±2 kpc 


(1) These are lower limits. 


Table 6 : Hi parameters ofNGC 253. 


erloo, Fraternali & Sancisi 2007) and NGC 2403 (Fratemali et al. 

2001 ). 

In order to isolate the anomalous Hi emission in NGC 253 
from the regular Hi, PV slices aligned with the major axis were 
extracted along the full extent of the galaxy’s minor axis. Each PV 
slice was visually inspected. Any emission deemed to be kinemati¬ 
cally anomalous was isolated, removed, and used to construct a new 
cube. In Figure 8 , the anomalous Hi is delimited by the white con¬ 


tours. This position-velocity diagram is representative of the ma¬ 
jority of PV slices with anomalous Hi. 

The spatial distribution of the anomalous Hi in that new cube 
is delimited by red contours in the channel maps shown in Fig¬ 
ure 9. Summing all the gas gives an extended component on the 
approaching half of the galaxy and a smaller component on the 
receding side as can be seen in Figure 10. The mass of that com¬ 
ponent is estimated to be ~3.5% of the total Hi mass of the galaxy 
(~ 7.8 X 10^ Mq). This is equivalent to the mass of the anomalous 
gas found by Boomsma et al. (2005). 

Naturally, this should be seen as a lower limit since our tech¬ 
nique has most probably missed a large fraction of the "anomalous" 
gas, especially along the minor axis to the NW and SE due to pro¬ 
jection effects. This can be seen on the PV slices along the minor 
axis shown in Figure 11. However, this will not be a problem for 
the kinematical study since, as we will see in Sec. 6 , the data along 
the minor axis will be excluded. 

In order to generate a velocity field representative of the circu¬ 
lar kinematics of the Hi disk, we again used PV slices aligned with 
the major axis of the galaxy, and fitted interactively three Gaussians 
to all of the line profiles making up a slice. Great care was taken in 
selecting the Gaussian associated with the disk of the galaxy. This 
was done by simultaneously viewing a position velocity slice and 
the fitted Gaussians. The line profile along the white dashed line on 
the left part of Figure 8 is shown on the right. For this particular 
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Figure 8. (left) PV slice through the centre of the galaxy. The velocities corresponding to the peaks of the Gaussians have been 
indicated with the horizontal dashed lines, while the vertical dashed white line shows the location of the profile on the right, (right) 
The red Gaussian traces the HI ridge, the orange Gaussian traces the beam-smeared emission, and the green Gaussian is fitting the 
anomalous HI. 


line profile, the red Gaussian parameterizes the radial velocity of 
the high-intensity HI ridge in the PV slice. This is the component 
of the line profile associated with the regularly-rotating Hi disk of 
the galaxy. The orange Gaussian, with its peak shifted toward sys¬ 
temic velocity, represents the beam-smeared component of the line 
profile. The green Gaussian captures the kinematically anomalous 
Hi emission. 

Stepping along the minor axis of the galaxy, this procedure 
was repeated for all the PV slices aligned parallel to the major axis. 
In doing so, a velocity field representing the circular kinematics of 
the HI disk of NGC 253 was built up pixel-by-pixel. When none 
of the Gaussians fitted to a line profile provided a reliable measure 
of the circular rotation, a blank was assigned to the corresponding 
position in the velocity field. We consider the velocity field shown 
in Figure 12 (top panel) to better represent the circular kinematics 
of NGC 253 than the traditional intensity-weighted-mean velocity 
field. From the cube built for the anomalous Hi, a velocity field 
was also obtained by a moment analysis and is shown in the bottom 
panel of Figure 12. 


6 Hi KINEMATICS OF THE DISK AND HALO GAS 

To derive the RCs, we used the implementation of the tilted ring 
model (Rogstad, Lockhart & Wright 1974) in the GIPSY (van der 
Hulst et al. 1992) task ROTCUR (Begeman 1989, 1987). In the 
tilted ring model, a set of concentric rings is used to describe the 
motion of the gas in the galaxy. The gas is assumed to be in cir¬ 
cular motion. Each ring is characterized by a set of 5 orientation 
parameters, namely: a rotation centre (xc, TJc), a systemic velocity 
Vsys, an inclination i, a Position Angle PA and by a rotation ve¬ 
locity Vc- Naturally, the rotation centre (xc, j/c) and the systemic 
velocity Vsys should be the same for all the rings, at least within 
the optical disk, but i and PA will vary if the Hi disk is warped 
outside D 25 . 

The line of sight velocity at any (x, y) position in a ring with 
radius R is given by 


V (x, y) = Vsys + Vcsin{i)cos{0) (2) 

where 9 is the position angle with respect to the receding major axis 


V (km s *) 
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Figure 12. Velocity fields of the disk (top) and anomalous (bottom) Hi. 


measured in the plane of the galaxy. 0 is related to the actual PA 
in the plane of the sky by 

-{x-xo)sin{PA)p{y-yo)cos{PA) 
cos(0) = - — - (3) 


-{x-xo)cos{PA) + {y-yo)cos{PA) 

^ ~ ' Rcos{i) ^ ' 
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R.A. (2000.0) R.A. (2000.0) R.A. (2000.0) 

Figure 9. Individual channel maps identifying the anomalous HI emission.The solid black contours are at levels of 0.8 0.6 0.5 0.4 0.3 
0.2 0.1 0.05 0.04 0.02 0.01 0.006 0.003 Jy/beam. The dashed black contours are at -0.1 -0.05 -0.01 -0.003 Jy/beam. The red contour is 
at 0.003 Jy/beam. The beam is shown in the bottom left corner of the first channel. 


6.1 Tilted-ring model for the Hi disk 

A \ cos9\^ weighting function and an exclusion angle of 60° about 
the minor axis have been used to give maximum weight to the ve¬ 
locity points close to the major axis and to minimize the influence 
of large deprojection errors close to the minor axis in view of the 
large inclination of the galaxy. We used two rings per beam size. 

The method consists of finding for each ring the best set of the 
5 orientation parameters (xc, Vc), Vsya, i and PA which minimizes 
the dispersion of Vc inside the ring. The following procedure is 
used: 


• The rotation center {xc,yc) and the systemic velocity Vays 
are determined first by keeping i and PA fixed (using the optical 
values). The rotation center and the systemic velocity have to be 
determined simultaneously since they are correlated. 

• Now, keeping (xc, Vc) and Vaya fixed, i and PA are fitted to 
map any possible warp of the Hi disk. The warps usually start just 
outside the optical. However, since the Hi disk in NGC 253 is about 
the same size as the optical, we may not reach the warp region. Here 
also, i and PA have to be determined simultaneously since they are 
correlated. 
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Figure 10. Combined WISE (RGB) + FUV (magenta) + X-rays (yellow: Rosat) + HI contours of the anomalous Hi (scale: l'~ 1 kpc). The contours correspond 
to 0.03, 0.2, 0.5, 1.0, 1.5 2.0 ly/beam km s"'^. 


• The previous two steps were done using the data of both sides 
of the galaxy together. Using the same fixed [xc, Vc) and Vsys, the 
previous step is repeated for the approaching and receding sides 
separately to identify possible departures from axisymmetry. 

The ROTCUR solutions (approaching, receding, both sides) 
are shown in Figure 13 with PA and i free to vary. Because of 
the low resolution of our data, it can be seen that, while PA is 
fairly well determined (< PA >= 235° ± 5°, i is not constrained 
at all. It is clear, looking at the galaxy, that inclinations ~ 40°, 
as seen for r > 10', are completely unrealistic. However, much 
higher resolution data (30”) exist from Fabry-Perot observations 
(Hlavacek-Larrondo et al. 2011). In their Figure 10, we can see 
that the inclination is nearly constant at i = 76° ±4°. Additionally, 
similar values of i (78.5°) and PA (232°) have been inferred from 
CO observations (e.g. Sofue et al. 1999; Hougton et al. 1997). In 
the case of NGC 253, adopting constant PA and i is reasonable 
since, as we have seen in Sec. 4, the Hi disk has the same size as 
the optical disk and warps usually start past the optical radius. We 
thus adopt constant values of PA = 235° and i = 76° (green lines 
in Fig. 13) to derive our adopted RC. 

A similar approach was used by Puche, Carignan & van 
Gorkom (1991) using higher resolution Hi data. The adopted RC 
for the Hi disk (solution for both sides with PA and i fixed) is 
given in Table 7. An important point to notice is the large errors 
on the first two points in the inner parts. This is expected for at 


least three reasons: 1- the strong absorption feature in the centre; 
2- the small quantity of Hi in the centre (cf. Fig, 7); 3- the strong 
disturbance expected from the starburst in the inner regions. 


6.1.1 Three-dimensional modelling 

As a further test on our choice of i, three different tilted ring models 
were generated, with constant inclinations of 66, 76 and 86 degrees. 
For each of the models, all the other parameters (PA, Vays,Xc, Yc) 
were kept fixed to the same constant values, allowing only 14 to 
vary with radius. The tilted ring models were then used to gener¬ 
ate 3D models of the NGC 253 Hi data cube. This was done using 
the task GALMOD in GIPSY. For the radial distribution of Hi we 
adopted the Hi surface density profile shown in Figure 7. A fiducial 
Hi disk scale height of 0.5 kpc was used for each model. The re¬ 
sulting model cubes were smoothed to the spatial resolution of the 
data, and the fluxes scaled based on the NGC 253’s integrated flux 
density of 728 Jy km/s. 

These models, each based on a different inclination, were then 
compared to the data on a channel-by-channel basis. Doing so 
clearly showed that the i — 76° model consistently matches the 
data better than the models with i = 66° and i = 86°. For three 
different channels. Figure 14 shows the comparisons between the 
models and the data. Clearly, i = 66° (left) produces emission too 
elongated in the minor axis direction and i = 86° (right) too nar- 
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Figure 11. PV slices parallel to the major axis. One is on the centre while the other two are 4' NW and SE of the major axis along the minor axis. Contour 
levels are -3, 3, 6, 12, 24, 48, 96 mJy/beam 


row. Based on these results, we adopt i = 76° as the preferred 
inclination for the Hi disk of NGC 253. This inclination, together 
with the rotation curve from the tilted ring model, accurately de¬ 
scribes the kinematics of the HI disk. 


6.1.2 Comparison of the RC with the data 

Finally, the best test to check that our RC accurately represents the 
kinematics of the disk is to look at the superposition of the reced¬ 
ing and the approaching sides on a Position-Velocity (PV) diagram 
obtained along the major axis, as shown in Figure 15 (top panel). It 
can be seen that the deprojected velocities follow very closely the 
ridge of the emission, which shows that our RC is a good represen¬ 
tation of the kinematics on the major axis. Note on the PV diagram 
the anomalous Hi around 200 km s“^. We should point out the 


presence of Galactic Hi, which can be seen around 0 km s“ ^. Even 
though we did exclude the channels with obvious Galactic emis¬ 
sion, it can be seen that some of it is still present. However, since 
there is much more Hi on the receding side than on the approaching 
side, very few data points from the approaching side contribute to 
the last points of the solution for both sides (cf. Fig. 13) and make 
us confident in our adopted RC. 

One last point needs to be emphasized. At the beginning of 
the section, we mentioned that a \cos9\‘^ weighting function and an 
exclusion angle of 60° about the minor axis have been used to min¬ 
imize the influence of large deprojection errors close to the minor 
axis in view of the large inclination of the galaxy. Typically, when a 
galaxy is not highly inclined (p.e. i = 30 — 45°), all the data can be 
used. For intermediate inclinations (p.e. i = 45 — 60°), a free angle 
of ~ 30° and a |cos0| weighting function should be used to give 
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Figure 13. ROTCUR solution and derived RC for the Hi disk.The green 
lines show our adopted PA = 235° and i = 76°. 


66 deg 76 deg 86 deg 



R.A. (2000.0) R.A. (2000.0) R.A. (2000.0) 


Figure 14. Comparison of the data (black) with 3D-models (red-dashed) 
using ROTCUR solutions with i = 66° (left), 76° (centre) and 86° (right) 
for 3 different channels. The contours are at levels of 0.03, 0.2, 0.4, 0.7 
Jy/beam. 
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Figure 15. PV diagrams along the major axis with the receding (red pluses) 
and approaching (blue pluses) RC superimposed using the adopted kinemat- 
ical parameters with a free angle of 60° and a \cos9\'^ weighing function at 
the top and no free angle and no weighing function at the bottom. The black 
contours are negative contours around the absorption feature. 


more weight to the data points close to the major axis. However, 
for high inclinations (p.e. i > 75°), it is suggested to use a larger 
exclusion angle (~ 60°) and a \cos6\^ weighting function. 

To illustrate this, we ran a ROTCUR model for NGC 253 with 
no exclusion angle and no weighting function. The resulting RCs 
for the approaching and receding sides are overlaid on the PV- 
diagram of the major axis in the bottom panel of Figure 15. We 
see that this solution is inconsistent with the kinematics observed 
on the major axis and must be discarded. This illustrates well that 
one has to be very careful when deriving the kinematics of highly 
inclined galaxies, especially with ROTCUR. 
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Radius 

(arcsec) 

V,at 

(kms“^) 

AV 

(km s“i) 

110 

79.3 

61.6 

220 

183.3 

33.1 

330 

198.3 

8.6 

440 

196.2 

8.0 

550 

203.7 

6.4 

660 

211.4 

8.7 

770 

214.1 

10.1 

880 

211.8 

13.2 

990 

198.6 

15.5 

1100 

186.7 

16.5 


Table 7: Rotation velocities and errors for the Hi disk of NGC 253. 


0 


Radius (kpc) 

8 _ TZ 16 


20 


300 : 

250 - 

7 200 - 

co 

a 

> 

101 


^ Disk HI 
^ Puche et al. (1991} 




♦ t 




0 - 


% ' ' ' 200 ' ' ' 400 ' ' 'eoo' ' 'soo' ' '1000 ' '1200 
Radius (arcsec) 


Figure 16. Comparison of the VLA & KAT-7 rotation curves. 
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Figure 17. ROTCUR solution and derived RC for the anomalous Hi. The 
green lines show our adopted PA = 230° and i = 76°. 


Figure 16 compares our adopted RC for the disk with the one 
derived by Puche, Carignan & van Gorkom (1991) using VLA ob¬ 
servations. We see that out to r ~ 12 kpc (the last measured velocity 
point with the VLA data), both rotation curves agree well within the 
errors, with the exception of our velocity at 5.5', which is ~ 15 km 
s“^ larger. As for the KAT-7 RC, it rises to a velocity of ~ 200 km 
s“^at 5.5', stays more or less flat out to ~ 15' and than declines 
by ~ 20 km s-^out to the last point. We will discuss this possible 
decline in Sec. 7.1. 

6.2 Tilted-ring model for the anomalous Hi 

The kinematical solution for the halo gas is shown in Figure 17. 
Tilted ring models were fitted to the velocity field of the anomalous 
HI (shown at the bottom of Fig. 12). When allowed to freely vary 
with radius, it is seen that PA and i are much better constrained 
for the anomalous component than for the disk. Despite the small 
change of PA for the last two points, it was decided to adopt con¬ 
stant PA and i (green lines). The rotation velocities for the anoma¬ 
lous gas are given in Table 8 and are shown in blue in Figure 18, 
along with the rotation velocities for the disk. It is different from 
the RC seen in the bottom panel of Fig 17 because PA and i were 
kept fixed (green lines). 

This RC for the anomalous gas is surprisingly very regular. 
It rises out to ~15 kpc and then declines slowly out to the last 
measured velocity point. In the outer parts, the anomalous gas is 



Figure 18. Comparison of the disk kinematics and of the anomalous Hi 
kinematics. 


kinematically lagging the disk gas by ~ 100 km s“^. Such lagging 
of extra-planar gas is also seen in other galaxies (see e.g. Elson 
2014; Hess et al. 2009; Heald et al. 2007; Fraternali et al. 2001). 
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Radius 

(arcsec) 

V,at 

(kms“^) 

AV 

(km s“i) 

440 

32.6 

16.9 

550 

55.7 

15.8 

660 

76.1 

15.4 

770 

96.3 

15.0 

880 

107.8 

18.5 

990 

98.6 

20.8 

1100 

91.9 

20.7 

1210 

80.7 

22.1 

1320 

75.2 

19.5 

1430 

70.2 

19.4 


Table 8 : Rotation velocities and errors for the anomalous Hi of 
NGC 253. 

7 DISCUSSION 

7.1 Is the Rotation Curve of NGC 253 really declining ? 

Since the discovery in the 1970s that the RC of most spiral galax¬ 
ies remains flat out to their last measured point (Freeman 1970; 
Bosma 1981), the obvious question is: where do the disks of spi¬ 
ral galaxies end? Especially in the case of spiral galaxies, the Hi 
surface density declines with increasing radius to a point where it 
decreases abruptly and the Hi disk seems to be truncated. This trun¬ 
cation is usually around surface densities ~10^®-10^® cm~^ (Hi 
self-shielding limit). This does not necessarily mean that there is 
no more hydrogen further out. As suggested hy Bland-Hawthom, 
Freeman & Quinn (1997), it may just be that the amhient radia¬ 
tion field out there is sufficiently strong that cold gas gets ionized 
(Maloney 1990; Silk & Sunyaev 1976). While the UV background 
radiation may not be sufficiently powerful to be responsible for this 
ionization, the situation may be different in a starburst object like 
NGC 253. 

We have seen in Sec. 4 that the Hi disk of NGC 253 only ex¬ 
tends out to the optical radius and the WISE mid-IR extent. But 
if, like most late-type spirals, there is hydrogen further out and if, 
also as seen in most spirals, the cold gas disk is warped, it could 
he exposed to the intense UV radiation of the starburst nucleus and 
be ionized. This was the motivation for Bland-Hawthorn, Freeman 
& Quinn (1997) and Hlavacek-Larrondo et al. (2011) to try to de¬ 
tect Ha emission in the outer parts of NGC 253 using Fabry-Perot 
(FP) interferometry, and perhaps succeed to find where the mass 
(luminous & dark) ends. 

Both groups detected Ha at the limit of the Hi disk (at 690, 
720 and 900 arcsec) but Hlavacek-Larrondo et al. (2011) also de¬ 
tected [Nil] out to 1140 arcsec or 1.4 D 25 on the receding side 
of the galaxy. The most interesting result was that both sets of 
measurements imply that the RC is declining past the last previ¬ 
ously measured Hi velocity point (Puche, Carignan & van Gorkom 
1991). From this, Bland-Hawthom, Freeman & Quinn (1997) con¬ 
clude that it suggests that the dark halo of NGC 253 may be tmn- 
cated near the Hi edge and provides further support for the link 
between dark matter and Hi (Bosma 1978; Carignan & Beaulieu 
1989; Carignan & Puche 1990a,b; Puche, Carignan & Bosma 1990; 
Jobin & Carignan 1990; Hoekstra, van Albada & Sancisi 2001; 
Swaters et al. 2012; Meurer, Zheng & de Blok 2013). Most impor¬ 
tantly, it would provide a means to probe the gravitational potential 
beyond the edge of the Hi disk. 

Figure 19 compares our Hi RC to the optical FP velocity mea- 



Figure 19. Compaiison of the Hi disk kinematics with the Fabry-Perot data 
of Bland-Hawthom, Freeman & Quinn (1997) and Hlavacek-Larrondo et 
al. (2011). 

surements on the receding side of NGC 253. Within the errors the 
Hi and optical velocities agree very well and seem to confirm that 
the RC of NGC 253 declines for R >12 kpc. Could a bad choice 
of orientation parameters, mainly PA and/or i, mimic a declining 
RC? As seen in Fig. 13 (see also Koribalski, Whiteoak & Houghton 
1995), there is a suggestion that PA might be increasing for the last 
two rings. Tilted-ring models, mn with those higher values, did not 
change significantly the derived velocities. However, i could have 
a more important effect but to bring the last point at the level of the 
flat part of the RC, an i ~ 56° would be needed which, as shown 
in Sec. 6.1.1, is excluded. 

NGC 253 is the second galaxy in the Sculptor group with a 
declining RC. Carignan & Puche (1990a) and de Blok et al (2008) 
showed NGC7793 to also have a declining RC, which is contrary 
to what is observed for the vast majority of spiral galaxies. 

7.2 Star formation in NGC 253 using WISE data 

If we want to understand the gas flows in the starburst galaxy NGC 
253, it is important to understand its star formation. WISE imag¬ 
ing provides the ideal platform for global property studies of large 
galaxies, such as NGC 253, due to the large field of view and the 
broad photometric bands that cover both the stellar (3.4 and 4.6 fim, 
or the W1 and W2 bands) and ISM components (12 and 22 /rm, or 
W3 and W4 bands). Reconstructions using a "drizzle" sampling 
technique (Jarrett et al. 2012 ) that conserves the native angular res¬ 
olution (6 arcsec in the 3.4 ^m band) were carried out for NGC 
253. 

The resulting mosaics fully cover the galaxy and its immedi¬ 
ate environment. However, a significant impediment to the 12 and 
22 fim imaging arises from the nuclear starburst itself, saturating 
the inner few pixels that comprise the unresolved nucleus. To re¬ 
cover this lost data, we apply a PSF-technique that was developed 
by T. Jarrett for saturation in Spitzer IRAC and MIPS imaging (see 
IRAC Instrument Handbook) and successfully deployed for study 
of the Circinus Galaxy (For, Koribalski & Jarrett 2012). The recti¬ 
fied WISE imaging of NGC 253 is shown in Figure 5 and Figure 
10, where all bands are combined into a four-color graphic: W1 
assigned shades of blue, W2 green, W3 orange and W4 red. 

Deploying the characterization pipeline of the WISE En¬ 
hanced Resolution Atlas (Jarrett et al. 2013), foreground stars are 
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identified and removed, the two-dimensional shape of NGC 253 
is determined and photometry is extracted. The one-o- isophotal 
(23.24 Vega mag arcsec“^ or 26.6 mag arcsec“^ in AB) diam¬ 
eter is found to be 42.2 arcmin with an axis ratio of 0.276. The 
corresponding integrated flux densities are 12.74, 8.44, 47.95, and 
118.04 Jy, for Wl, W2, W3 and W4 respectively, indicative of very 
strong PAH emission (W3) and warm dust continuum (W4) arising 
from the active star formation. We note that the "total" flux, as in¬ 
ferred from the isophotal plus extrapolated disk emission, is only 
a few per cent larger, consistent with a relatively truncated stellar 
disk. 

As demonstrated in Jarrett et al. (2013) and Cluver et al. 
(2014), the stellar mass-to-light ratio has a simple dependence on 
the WISE Wl and W2 integrated fluxes. With our adopted distance 
of 3.47 Mpc, the corresponding Wl absolute magnitude and Log 
[in-band luminosity L^i (L©)] is -24.24 and 10.99 , respectively. 
The W1-W2 color is 0.21 mag, and employing the M/L relation 
from Cluver et al. (2014), the inferred Log [stellar mass (Mq)] is 
then 10.33, indicative of the significant (mass) evolved stellar pop¬ 
ulation that represents the "back bone" of NGC 253. 

Star formation activity may be gauged by the WISE 12 and 
22 ^m luminosities. Log [L (L©)] = 9.67 and 9.77, respectively. 
Here we apply the relations in Cluver et al. (2014) to estimate the 
dust-obscured star formation rate (SFR): 4.9 and 5.1 Mq yr“^, re¬ 
spectively based on the 11.3 ^m PAH emission and the 22 /rm dust 
continuum. Although this global SFR is typically small compared 
to luminous infrared galaxies, it is significant in the Local Volume, 
similar to the starburst M82 and to the barred grand-design spiral 
M83 (Jarrett et al. 2013). 

We infer a star formation density of 1.7 Mdot/year/kpc^ for 
the central kpc of NGC 253, consistent with a nuclear starburst. For 
normal galaxies, Cluver et al. (2010) showed a clear linear trend in 
global SFR relative to the Hi gas content, effectively the Kennicutt- 
Schmidt law. Starburst galaxies, such as M82 and NGC 253, devi¬ 
ate from this relation due to the enhanced nuclear starburst phase, 
as seen in Figure 20. Finally, normalizing the SFR with the stellar 
mass, the resulting Log [specific SFR (yr~^)] is -9.86 and -9.83, 
respectively W3 and W4. This disk building rate is relatively high 
for nearby galaxies - NGC 253 is still actively building its disk 
and bulge D and is comparable to other large star forming spirals, 
notably M83 and NGC 6946 (Jarrett et al. 2013). 

7.3 Origin of the Anomalous Hi: inflow or outflow ? 

As discussed in the introduction, at least three possible scenarios 
could be envisaged to explain the anomalous (halo) gas that is ob¬ 
served. The first one is the classical galactic fountain outflow sce¬ 
nario, where gas is being expelled from the disk through multiple 
supernova explosions from clusters of massive stars (Schapiro & 
Field 1976). The second one is the inflow scenario of gas of exter¬ 
nal origin in the form of primordial gas clouds left over from the 
formation of the galaxy (Oort 1966). Finally, that halo gas could 
come from ISM tom out of dwarf galaxies during a close encounter 
with NGC 253 (Putman, Peek & Joung 2012). 

It is clear from this study that the first scenario can explain 
a large part of the anomalous gas observed. The known starburst 
nature of NGC 253, with its SFR > 5 M© yr“^ (see Sec.7.2) and 
the large amount of hot halo gas seen along the minor axis in the 
X-ray observations (Pietsch et al. 2000), all point to gas being ex¬ 
pelled from the nuclear region. This starburst-driven super wind 
was studied in great detail by Westmoquette, Smith & Gallagher 
(2011) through deep Ha imaging with the MPG/ESO 2.2 m WFI 



Figure 20. Obscured star formation rate compared to the total neutral hy¬ 
drogen content. Adapted from Cluver et al. (2010), the sample is from 
SINGS (Kennicutt et al 2003) with highlighted galaxies, notably nearby 
starbursts M82 and NGC 253 (this work). 

and optical spectroscopy with the VLTWIMOS-IFU and with the 
WIYN/SparsePak IFU. They investigate the known minor axis out¬ 
flow cone, which is well-defined in the Ha imaging and kinemat¬ 
ics between radii of 280 and 660 pc from the nucleus. Kinematic 
modelling indicates a wide opening angle (~60°), an inclination 
consistent with that of the disc and deprojected outflow speeds of a 
few 100 km s“^ that increase with distance above the plane. 

A dual origin for the gas in the halo in the form of galac¬ 
tic fountains originating within the boiling star-forming disk and 
a strong galactic super wind emanating from the nuclear starburst, 
has previously been suggested by Pietsch et al. (2000). It is thus 
quite possible that the classical galactic fountain (ionized outflow 
seen in the Ha emission, cooling down through radiative processes 
and raining back onto the disk) provides most of the extra planar 
gas up to 5-10 kpc, while the gas higher up in the halo could origi¬ 
nates from the starburst region, expelled as hot gas and then cooled 
to its neutral form.The most convincing piece of evidence that sup¬ 
ports the outflow origin for the extra planar gas is seen in Figure 
12 , where it is apparent that the halo gas retained the kinematical 
signature of the disk. 

As suggested by a referee, the location of the anomalous gas, 
situated more toward the edge of the disk (see Fig. 10) could pos¬ 
sibly be interpreted as the result of the superwind outflow, as in the 
case of NGC 1482 (Hota & Saikia 2005). It looks as if the anoma¬ 
lous gas in the halo is pushed by this outflow and is providing colli- 
mation to the outflowing gas. This halo gas would have a different 
origin than the outflow gas, most probably through galactic foun¬ 
tains, which is suggested by the observation that it retains traces of 
the main disk kinematics. 

If the anomalous gas was an inflow from external origin, it is 
expected that those clouds would include material at forbidden ve¬ 
locities. After a careful visual search of the Hi cube of NGC 253, no 
such clouds with that expected peculiar kinematics were uncovered. 
However, looking at the size of the inflow clouds, e.g. in NGC 891 
(Oosterloo, Fraternali & Sancisi 2007), of only 1-2', such clouds 
would most probably not be detected in our data cube because of 
beam dilution (beam ~ 3.5' x 3'). Thus, we cannot exclude com- 
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pletely the presence of such small inflow clouds. Finally, our data 
seems also to exclude gas torn off a dwarf galaxy during a close 
encounter since again we would not necessarily expect that gas to 
have the same kinematics as that of the disk. 

The last question left to answer: where spatially is most of the 
detected anomalous gas? Looking at Fig. 18, it can he seen that the 
halo gas has a much shallower gradient than the disk gas and is 
lagging by ~ 100 km s“^ in the outer parts. Modelling done (see 
e.g. Fraternali et al. 2005) of extra planar gas predicts a lag of ~15 
km s“^ kpc“^, consistent with observations of edge-on galaxies. 
Applied to NGC 253, this suggests that the bulk of the gas should 
be around 6-8 kpc from the plane of the main Fll disk. If it were 
closer to the disk (see e.g. Hess et al. 2009), it would have a steeper 
gradient in the inner parts and a much smaller lagging velocity. 


7.4 Origin of the starburst nature of NGC 253 

One of the principal reasons often given for high SFR is galaxy- 
galaxy interaction (see e.g. Larson & Tinsley 1978), whether in 
the form of merger, close encounter or interaction with the IGM. 
A good example is the M81 group where the three core members 
M81, M82 and NGC 3077 are known to be closely interacting from 
Hi studies (van der Hulst 1979; Yun, Ho & Lo 1994) and, at the 
same time, all three show evidence for an AGN or starburst activity, 
most likely induced by the on-going tidal interactions. 

In contrast, the environment of NGC 253 appears much more 
devoid of any sign of gravitational interaction either with low-mass 
companions or the IGM. The search of Cote et al. (1991) for dwarf 
galaxies in the Sculptor group (see their Fig. 4) do not show any 
obvious nearby dirr with Hi emission, while a few possible dSphs 
candidates were identified. As for the presence of a significant IGM 
component, the study by Putman et al. (2003) suggests that the sur¬ 
rounding region has very little intergalactic gas (see their Fig. 12). 

However, the apparent lack of evidence of any actual interac¬ 
tion does not tell us anything about possible past interactions of 
NGC 253 with its environment. For example, as shown by Whiting 
(1999), the spin of the NGC 253 disk is consistent with torquing 
from NGC 247, which is its closest (~ 350 kpc) large companion. 
Moreover, both galaxies appear to have truncated Hi disks (Puche, 
Carignan & van Gorkom 1991; Carignan & Puche 1990b), since 
they barely extend out to their optical diameter. This is contrary 
to what is generally seen in late-type spirals where the Hi disk is 
usually much more extended than their stellar disk (Sancisi 1983). 
As an example, the Hi disk of NGC 300, another late-type spiral in 
the Sculptor group, has a diameter more than 1.5 times the optical 
diameter (Puche, Carignan & Bosma 1990). 

Despite its apparent isolation, there is some evidence which 
suggests that NGC 253 may have been recently involved in a 
merger. For example, in a recent X-ray study to spatially resolve 
the starburst region of NGC 253 (Wik et al. 2014), the 4-6 keV data 
clearly shows what appears to be a double nucleus (see their Fig. 
4). Earlier studies (Anantharamaiah & Goss 1996; Prada, Gutierrez 
& McKeith 1998) also found evidence of two dynamically distinct 
systems close to the centre of NGC 253. As discussed by Davidge 
(2010), one system has an axis of rotation that differs from that 
of the galaxy disk, while the other appears to be counter-rotating 
with respect to the galaxy disk. These are all signatures of a re¬ 
cent merger. Das, Anantharamaiah & Yun (2001) estimate that the 
merger may have involved the accretion of 10® M© of material 
~10^ years ago. In view of the actual isolation of NGC 253, this is 
the most likely explanation for the starburst activity in its nucleus. 


8 SUMMARY AND CONCLUSIONS 

This study presents the results of ~150 hours of KAT-7 observa¬ 
tions in Hi line mode, of which ~115 hours were spent on NGC 
253. For the analysis, the data were smoothed to a velocity reso¬ 
lution of 5 km s“^ and the spatial resolution (synthesized beam) 
is 213"x 188". As much as 33% more flux was detected by KAT- 
7, compared to previous VLA observations for a total Hi mass of 
2.1 ± 0.1 xlO® M0. A sensitivity limit of ~L0 x 10^® cm“® in 
column density allows us to uncover a large quantity of extra planar 
gas out to projected distances from the plane of ~9-10 kpc in the 
centre and 13-14 kpc at the edge of the disk. However, in the radial 
direction, the Hi disk is relatively small for a late-type spiral, sim¬ 
ilar to the optical and mid-IR disk, consistent with the tmncation 
scenario. 

A robust method was developed to separate the regular and 
anomalous Hi components in NGC 253. This method, which used 
interactive profile fitting in combination with PV diagrams, allowed 
to fit separately the different components pixel by pixel. For the 
disk, the RC derived from the KAT-7 data agrees well with the VLA 
data out to ~12 kpc (the last measured VLA velocity point). For R 
>12 kpc, the RC declines by ~ 20 km s“^ out to the last point at 
~18 kpc. This decline, observed in the outer parts, agrees with pre¬ 
vious optical FP observations. As for the halo gas, which extends 
out to ~24 kpc, its RC has a very shallow gradient in the inner 
parts and is kinematically lagging the disk gas by ~100 km s“^in 
the outer parts. That extra planar cold gas component is seen at the 
edge of the hot gas component revealed by the X-ray observations. 

Our results point towards an outflow scenario for the origin of 
the extra planar gas. This outflow has two origins. First, as clearly 
seen in X-ray observations, a lot of hot gas is expelled in the halo 
from the nuclear starburst region. As the gas expands, it eventu¬ 
ally converts back to Hi as it cools through radiative losses; raining 
back down on the disk and feeding subsequent star formation. Sec¬ 
ondly, another part of that extra planar gas has most likely a galactic 
fountain origin from the star-forming disk. 

The KAT-7 observations presented in this work clearly show 
that despite its relatively small size (7 x 12 m antennae), this tele¬ 
scope occupies a niche for detecting large scale low surface bright¬ 
ness diffuse emission over the ~1° FWHM of its antennae. It 
should be kept in mind that this telescope was built primarily as 
a testbed for MeerKAT and the SKA such that any scientific result 
that can be obtained is a bonus. While most of the extragalactic Hi 
sources would be unresolved by the ~4' synthesized beam, many 
projects such as this one and the project on NGC 3109 (Carignan 
et al. 2013) can be done on nearby very extended objects such as 
Local Group galaxies or galaxies in nearby groups like Sculptor. 
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